All-sky incoherent search for periodic signals with Explorer 2005 

data 
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The data collected during 2005 by the resonant bar Explorer are divided into 
segments and incoherently summed in order to perform an all-sky search for periodic 
gravitational wave signals. 

The parameter space of the search spanned about 40Hz in frequency, over 23927 
positions in the sky. Neither source orbital corrections nor spindown parameters 
have been included, with the result that the search was sensible to isolated neutron 
stars with a frequency drift less than 6 • 10 _11 Hz/s. 

No gravitational wave candidates have been found by means of the present analysis, 
which led to a best upper limit of 3.1 • 10 -23 for the dimensionless strain amplitude. 
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PACS numbers: 95.55Ym, 04.80.Nn, 95.75.Pq, 97.60. Gb 



A. Introduction 

The search for periodic gravitational wave signals is a stimulating challenge for data 
analists because of the considerable amount of computing time required. 

For blind searches, i.e. without any a priori knowledge about the source, a fully coherent 
analysis can not handle more than a few days of data because of the steep dependence of 
the size of the parameter space on the frequency resolution. 

In [l| three data sets, each two days long, from the Explorer 1991 run have been coherently 
studied by means of the F statistics method which led to and an upper limit of 1 • 1CT 22 
on h in the narrow band 921. 00-921. 76Hz. 

A similar technique, applied in [3] to the most sensitive 10 hours of the the LIGO S2 run, 
led to an upper limit of 6.6 • 10 -23 for isolated neutron stars in the band between 160 and 
728.8Hz. 

With the widening of the frequency band, due to the advent of interferometers as well as 
to improvements in the readout of resonant detectors, several incoherent and semi-coherent 
methods have been conceived and employed. 

n 

In (4j the Hough transform technique has been applied to the LIGO S2 data to perform 
a blind search for isolated neutron stars on a set of narrow frequency bands in the range 
200-400Hz, and a best upper limit of 4.43 • 10~ 23 has been set. 

In the present work, the simple technique of adding power spectra has been applied to 
the most sensitive 40Hz band of the 2005 run of the Explorer bar [51], resulting in a further 
improvement on the best upper limit on h, which is set to 3.1 • 10~ 23 at 920.14Hz. 

According to the results reported at the recent Amaldi7 conference, the analysis of the 
LIGO S4 run 6j is leading to a sensible improvement in this direction (about an order of 
magnitude); remarkably, the limit set in the present work is still competitive with the one 
coming from the S4 data in the same frequency band. 
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Figure 1: The noise level of Explorer during the 2005 run. The solid horizontal line is the cut 
applied to select the best spectra. 

B. The data set 

At the end of April 2005, after a short commissioning break, the resonant antenna Ex- 
plorer was again on air and operated with the usual, remarkable duty cycle (86% from April 
to December 2005) and a good stability. 

The data stream taken by the bar until the end of 2005 (after which the sensitivity curve 
has been modified) has been divided into 25161 segments, each about 14 minutes long. 
The most sensitive band of the Fourier transform of these segments, namely Nf = 32178 
frequency bins in the range 885-925Hz, has been selected for the analysis. 

A noise cut has been applied to the total power contained in each spectrum, with the 
purpose of discarding the noisy ones thus creating an homogeneous set of spectra. This 
allowed us to apply for the subsequent analysis the simple power addition method, without 
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Figure 2: The square root of the average spectrum represents the typical Explorer sensitivity curve 
during the 2005 run. 

weighting each spectrum with the corresponding noise level. 

This selection led to the creation of a data set {Si} made of the Ni = 11749 cleanest 
spectra (corresponding to an effective data time of 114 days), and of a second set containing 
just the best N2 = 3875 (used to deal with the critical zones of the spectrum, corresponding 
to the resonant modes of the bar, around 888Hz and 920Hz). 

The average sensitivity of the second data set is shown in Fig. [2] (the first set is similar 
except around the resonant modes); a few noisy lines may be noted, including small 1Hz 
harmonics on the left part. 



C. The analysis method 

For a given direction fj in the sky, the selected spectra have been deformed according to 
the Doppler shift formula 

/r * /exp (l - ^LJji) 
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and then summed and renormalized dividing by N\ or by N 2 : 



N ll2 



i=l 



The speed of the detector relative to the Solar System Baricenter has been computed 
thanks to the JPL ephemerides As the speed of the source have not been taken into 
account, the search is sensitive to isolated neutron stars, and not to those which are part of 
binary systems. 

Spindown has also been neglected: given the frequency resolution of 1.2 ■ 10 _3 Hz, this 
means being sensitive to an average frequency drift up to 6-10 _11 Hz/s during the observation 
period. 

The procedure has been repeated for any point of an optimized sky grid made of iV s k y = 23927 
possible directions, thus leading to the creation of a set {Sj} containing iV s ky "deformed and 
summed" spectra. 

The variance of the noise is obtained calculating, for each value of the frequency, the 
variance of the distribution of the iV s ky deformed spectra. The result, whose square root is 
shown in Fig. [3j agrees with the general expectation, based on the central limit theorem, 
that 



As expected, the plot shows an anomalous behavior of the noise variance in correspon- 
dence of the disturbances of the initial data set. These anomalous zones have not been 
taken into account for the candidate search, but have been included in the upper limit 
determination. 

The detection threshold is fixed by the requirement that the false alarm rate should be 
less than 1%. According to Poisson statistics one has to impose 



Si, 



(1) 




P(0, A) = e" A > .99, 



where the expected number of threshold crossings in absence of signal is 



A = p ■ N f ■ N sky , 
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Figure 3: Noise variance of the set of deformed and summed spectra. To help comparison with 
Fig. [21 the square root of a is actually shown. 

being p the probability of false detection in a single frequency bin and for a single direction 
in the sky, and Nf ■ N s k y the trial factor. 

One thus finds the condition p < 1.3 • 10~ n which, assuming that the iV s k y values of 
the shifted spectra at a given frequency bin are gaussian distributed, translates to a 7a 
threshold. 

In other words, a detection is claimed if, for some value of the frequency /, a given 
deformed spectrum Sj satisfies the condition 



Sj(f) - S(f) > 7a(f) 



(2) 



being S the average of the deformed spectra {Sj} as j spans over the N^y sky grid points. 

Since h oc y/S, to translate the detection threshold in h units one has to multiply the 
values shown in Fig. El by the factor 

_ '15 




where T is the length of each data segment and the factor 4/15 (the average angular sensi- 
tivity of the bar over the solid angle) is introduced to compensate the fact that amplitude 
modulation has not been taken into account when summing the deformed spectra. 
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Figure 4: Detection threshold (red, upper curve) and maximal triggers (blue, lower). The 1Hz 
disturbances on the left have been left only for illustrative purpose and cannot be considered as 
genuine threshold crossings. 

D. Results 



1. Cadidate search 

The threshold is shown in red in Figure HI while the blue line is the maximum value of 
h found, at any given frequency bin, among the set of the iV s ky spectra according to the 
following formula: 



^max (y ) 



max,- 



15 £,(/)-$(/) 



The anomalous regions of the spectrum have been cut out, since it is not reasonable to 
assume that they follow a gaussian distribution. To illustrate the point, the 1Hz disturbances 
have been keft in Figure HI to show that they would have produced fake candidates, if they 
had been included in the analysis. 

With these specifications, the maxima are never above the threshold, and thus no candi- 
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Figure 5: On the left, the outcome of the 16 injections, with initial amplitude given by the dashed 
red line. On the right, the upper limt curve. 



dates have been found. 



2. Software injections and upper limit 

To find an upper limit on h, a variation of the loudest event method 8| has been applied. 
This method allows to determine an upper limit starting from the loudest event present in 
a data stream, irrespectively of the fact that such an event may be due to noise or to a 
real signal. The idea is basically that if a strong real signal would have been present during 
the data taking, it would have produced an event louder than the loudest event actually 
recorded. This idea can be made quantitative by studying the detection efficiency of the 
experiment, for example by performing software injections at various SNR's. 

In our case, we have a loudest event for any frequency bin, and all these events form 
precisely the curve of maxima depicted in Figure H As a preliminary step, we had to 
determine our detection efficiency by injecting fake periodic signals in the Explorer data 
stream and finding out how hey did loook like after the analysis chain: the left plot of 
Figure [5] shows for instance the result of 16 injections of signals with h = 3 ■ 10~ 23 , equally 
spaced in frequency by 0.1Hz starting from 919.4Hz, and coming from randomly chosen 
directions in the sky. Then, the upper limit at 95% c.l. at has been determined for each 
frequency bin as the lowest injected amplitude which had produced a signal larger than the 
actual maximum at least in 95% of the cases. 

The right plot on Figure [5] shows the result, i.e. the curve of h upper limit at 95% 
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confidence level: the minimum is 3.1 • 10 23 at 920.14 Hz. 
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